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ABSTRACT 

The relativistic shock layer problem was numerically analyzed by using two relativistic Boltzmann- 
kinetic equations. One is Marie model, and the other is Anderson- Witting model. As with Marie 
model, the temperature of the gain term was determined from its relation with the dynamic pres- 
sure in the framework of 14-moments theory. From numerical results of the relativistic shock layer 
problem, behaviors of projected moments in the nonequilibrium region were clarified. Profiles of the 
heat flux given by Marie model and Anderson- Witting model were quite adverse to the profile of 
the heat flux approximated by Navier-Stokes-Fourier law. On the other hand, profiles of the heat 
flux given by Marie model and Anderson- Witting model were similar to the profile approximated by 
Navier-Stokes-Fourier law. Additionally we discuss the differences between Anderson- Witting model 
and Marie model by focusing on the fact that the relaxational rate of the distribution function depends 
on both flow velocity and molecular velocity for Anderson- Witting model, while it depends only on 
the molecular velocity for Marie model. 



I. INTRODUCTION 



Advances in elementary particle physics has made the study of the flow field of relativistic particles 
an important issue. The study of relativistic Boltzmann equations is particularly important for un- 
derstanding the fundamental properties of nonequilibrium relativistic gases. Despite this importance, 
numerical analysis of relativistic Boltzmann equations [TJ has not yet been reported. Fortunately, two 
relativistic kinetic equations can be analyzed numerically due to their simplified collision kernels. One 
is the Anderson- Witting model [2], and the other is the Marie model [3]. Both models are written in 
BGK (Bhatnagar-Gross-Krook) [Jj form. Assumptions in the construction of the numerical scheme 
to solve either the Anderson- Witting or Marie models include an inertial frame and no degeneracy 
of the particles. Projected moments [5] are calculated using a standard method, namely, the Eckart 
decomposition [BJ. 

In our analysis of the Marie model, we focused on the dynamic pressure, which is a characteristic 
quantity for relativistic gases. In the Marie model, the dynamic pressure is derived from the difference 
in temperatures of the gain and loss terms. The Anderson- Witting model avoids this inequality by 
using orthogonality of the four velocities of the flow by Landau-Lifshitz [8] to the nonequilibrium 
terms. However a concrete formulation of the temperature of the gain term has not yet been included 
in the Marie model. In our analysis, we therefore related the temperature of the gain term with the 
dynamic pressure in the framework of the N.S.F. (Navier-Stokes- Fourier) law derived from 14-moment 
theory (Tj [7] for the Marie model. Provided that we equate the temperatures of the gain and the loss 
terms in the Marie model, the bulk viscosity always becomes negative. 

To quantitatively examine the behavior of the Anderson- Witting and Marie models, we used a nu- 
merical method to solve the relativistic, steady-state shock-layer problem in the absolute standard of 
rest The meaning of "absolute standard of rest is not clear" . Please explain what you mean by this phrase., 
which is a hypothetical inertial frame. The heat flux and the dynamic pressure were determined from 
the simulations. The shock layer problem is suitable to investigate nonequilibrium gas dynamics in 
both the shock structure and in the boundary layer. 

Finally, we discuss differences between the Anderson- Witting and Marie models by comparing their 
flow and molecular-velocity dependent relaxation of the distribution function Do you mean "velocity 
distribution function" ?. Significant differences in relaxation of their distribution functions is observed in 
the negative tail structures The meaning of "negative tail structure" is not clear. Please explain what you 
mean, of their distribution functions at the shock's rising edge. Throughout this paper, the molecular 



potential is represented as a hard sphere molecule. 



This paper is organized as follows. In Section II we review the relativistic Boltzmann equation 
and its relation to the projected moments with Eckart decomposition. In Section III we describe the 
Anderson- Witting and Marie models, and describe a derivation of the temperature of the gain term 
from the N.S.F. law. In Section IV we describe a numerical method for solving the Anderson- Witting 
and Marie models. In Section V we describe an application of this method to solve the shock layer 
problem. 



II. RELATIVISTIC BOLTZMANN EQUATION 



The relativistic Boltzmann equation is written as [T] 

^ = w./> 

(fj'-f,f)Fadn^, (1) 
P*o 

where x a represents the four-dimensional coordinates, p a is a four-dimensional momentum vector, / is 
a distribution function defined by / = / ft, x l ,p l ) (i = 1, 2, 3), and F is the Lorentz invariant flux. In 
eq. (1), terms with prime indicate conditions after collisions and is the momentum space stretched 
by {^ 3 |(— oo, — oo, — oo) < (p l ,p 2 ,p 3 ) < (oo, oo, oo)}. x a , p a and F are given by 

x a = (ct,^ 1 ,^ 2 ,^ 3 ), (2) 

p a = mrf(v) (c, v , v 2 , v 3 ) , (3) 

p pQ p Q pl f ~ 2 1 ^ 

F = 50 = V(v-v*) vxv* 

c c V c z 



{plPa? ~ ™ 4 C 4 . (4) 



In eq. (3), 7(1)) is the Lorentz factor, which is given by 7(f) = — v 2 /c 2 . c is the speed of light 

and v l is the ith component of the particle velocity for i (= 1,2,3). In eqs. (3) and (4), m is the 
molecular mass. In eq. (4), g is M0ller's relative velocity. In eq. (1), a is the differential cross section 
and 0, is the solid angle. In eqs. (1) and (4), terms with an asterisk subscript belong to the collision 
partner. Rewriting eq. (1) in Lorentz variant form yields 

<>f + ^14 = / W - /*/) 9^dUdW (5) 



dt dx l 

Multiplying both sides of eq. (1) by p a and p a pP and integrating in momentum space, we obtain 



conservation equations in terms of N a = f^ 3 cp a f^- and T a/3 = J^ 3 cp a p^ f^Jr as 
Mass conservation: 

d a N a 

f a df d 3 p 

= C p 



m z dx a p° 
~ I (c + c-c-c)ff*Fadn^^ 

2J«3 xf 3 p [ i p U 

0. (6) 



Momentum-Energy conservation: 

d n T a/3 



a j3 d f d3 P 



' 1 n P ^dx a p° 



2 y^3 X ^ 3 v j pi p {> 

0. (7) 



According to Eckart [6], A*"" and T a @ can be decomposed to 



N a = c [ p a< ^- = nU a , (8) 

rraB f a B^P 

1 = CI P P P T- 

= p <a(3> - {p + w) A Q/3 



where n is the number density, 

p <af3> ig 

the shear stress, p is the isotropic pressure, w is the dynamic 
pressure, q a is the heat flux, e is the energy density, and U a is the four-dimensional velocity field of 
the flow given by 

u a = 1 (u){ c y), (io) 

where u l is the ith component of the flow velocity. A a @ in eq. (9) is the projector defined by 



where rj a " is given by 



v aP = nap 



Projected moments are obtained as [6] 
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Projected moments re, u l ,p <a ^ > , q a and w can be reduced from eqs. (14) and (16) to 14 projected 
moments n, u l ,p l;1 , q 1 , vj, = 1,2,3) as 



p<^> Ua = 0, 
q a U a = 0, 



(18) 
(19) 



Conservative equations in eqs. (6) and (7) yield balance equations for n (mass), U a (momentum), 
and e (energy) [T] as 

Dn + nV a U a = 0, (20) 



c 

4 ( /' ' 



DU p -wDU a -Dq a -q a VpU p - q p V p U a - ^U a q p DUp - U^^V pU^j 

(21) 

nDe = - (p + vj) V a U a + p <a ^> VpU a - V a q a + ^q a DU a (22) 



where D, V™, and the enthalpy per particle He are defined by 



D = U a d a , (23) 

V a = L a/3 - ^U a U^j dp = A a/3 dp, (24) 

h E = e + ^, (25) 
re 



where d a = -jj^a- Chapmann-Enskog expansion indicates that zu, p <a ^> and q a are approximated by 
the product of temporal-spatial gradients of projected moments and transport coefficients, the bulk 
viscosity r/, the viscosity coefficient u, and the thermal conductivity A as follows [TJ, 

w = - v V a U a , (26) 



P <af3> = 2/i 



\ (A-Aj + A? A?) - -A^ 



v 7 c/ 5 



(27) 

g a = A f V a # - -?-vV) (28) 
where 6 is the gas temperature. 

III. RELATIVISTIC KINETIC EQUATION 

In this section, we describe two kinetic equations: the Anderson- Witting and Marie models. 

Anderson- Witting model 

The Anderson- Witting model is given by [2] 

df _ U£p a ( f (0) 



^ = l/ w - /) . (») 

where /^°^ is an equilibrium function called the Maxwell-Jiittner function, and can be defined as 

^''■"'- wd^fe) '"^- (30) 

2 

where Qe is given by Qe = fe is the Boltzmann constant, and Oe is the temperature used in the 
equilibrium function f(°K K n is the nth order modified Bessel function and U£ is the four-dimensional 
velocity of the flow defined by Landau-Lifshitz [8] and written as 

U? = U a + (31) 
ne + p 

t in eqs. (29) is defined for the hard sphere molecule as pQ 
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(33) 



where v s is the relativistic speed of sound and G = K^(Q) / K2{Q- 

For the conservation law for N a and T a @ denoted by eqs. (6) and (7), the following constraints must 
be satisfied. 



(34) 



c / p a U La f^- = N a U La = c [ p a U La fW^ = N%U La , 
J3l 3 P Jffi P 

c [ p a p P U La f^ = TfUu, = c I pVM (0) ^ = Tfu La , (35) 

Jgi 3 P JM 3 P 

where subscript E indicates quantities derived from the equilibrium distribution function /(°). These 
constraints are considered to be satisfied by the orthogonality of U" to nonequilibrium terms in either 
N a or T a ". Multiplying eq. (35) by U Lf3 /(nc 2 ), we obtain [T][S] 

e E = e. (36) 

This relation in eq. (36) originates from the relation e = mc 2 (g(Q — j-J , 

9e = 0- (37) 

In the Anderson- Witting model Be is the temperature used in the equilibrium function f(°\ which is 
equal to 8, the temperature of the gas. 

Marie model 



The Marie model [3] is obtained by replacing L $ a in eq. (29) by m. The Marie model can be 
written as 

where r is defined in eq. (32). 

By multiplying both sides of eq. (38) by p^p 1 and integrating in momentum space, we obtain 

8 a T a ^ = — (T^ - T&) , (39) 

.Qj3„7 f d3 P 



where T Q ^ 7 = J^ 3 p a p^p y f —tr and can be decomposed to yield [T] 



r 2 / 
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+C A [p <af5> W + p <a ^>U p + p < ^ > [/°) . (40) 
C\ , C2 , C3 and C4 are functions of C shown in pQ . 

Multiplying both sides of eq. (39) by UpU y It is unclear what you are trying to say with "with eqs. (17) 
and (40)". and eliminating terms with nonequilibrium projected moments, we obtain 

e E -e = -^W, (41) 
n 

where ip is defined in Appendix A. From eq. (41), the energy per particle e is not conserved in the 
collision term in the Marie model. 

Multiplying both sides of eq. (39) by Ap^ with eqs. (17) and (40) and eliminating terms with 
nonequilibrium projected moments, we obtain 

™ = -V«E)V a U a = - (fj(( E ) + 7?(Ce)) V a C/ a . (42) 

In eq. (42), — J7(Cs)V a u a is the dynamic pressure derived from either pe~P or eg — e and —fj{C,E)^ a U a 
is the dynamic pressure derived from the left hand side of eq. (39). If ( = (e, — fj((E)V a u a = 0, 
w = —fj(()V a U a . 77(C) and t]((e) are given in Appendix A and are plotted in Fig. 1. Figure 1 shows 
that ry(C) is negative for all ranges of C (0 < C)- 
From eqs. (40) and (41), eg is given by 

eE = e H t-tEJ. (43) 

From the following approximate relation between e and p [T] , 

PF -p = n(e E -e) nk(e E - e) 



we obtain Qe or 9e as 



j_ _ 1 W(Ce) u ,, 

Ce C nmc2C v (( E )v((E) ' 1 ' 

nC v {e E )ri{pE) 



where C v is the const ant- volume specific heat. 



IV. NUMERICAL METHOD 



In this paper, we exclude photons, whose mass is zero, and molecules with velocities of the speed 



of light. For m/0 and v < c, the Anderson- Witting model defined in eq. (29) can be rewritten as 



ne + p 



(c, -«•) 



i\T 



(/ (0) " /) 



c 2 r 



(47) 



In general, the distribution function is / = f(t, x 1 , x 2 , x 3 , v 1 , v 2 , v 3 ), which has a one-to-one corre- 
spondence to / = f(t, x 1 , x 2 , x 3 ,p l ,p 2 ,p 3 ). In this work we use / = f(t, x 1 , x 2 , x 3 , v 1 , v 2 , v 3 ) instead 
of / = /(t, x 1 , x 2 , x 3 ,^ 1 ,^ 2 ,]? 3 ). This transformation is shown in eq. (47) and can be done readily for 
the Marie model. To calculate the projected moments, we transform d 3 p/p° into velocity space d 3 v 
as 



d 3 p 
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dp 1 
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m 2 j(v) 4 3 



d 6 v 



(48) 



From eq. (48), N a can be rewritten as 



N° 



c [ p a fd 3 p/p° 



c I mj(v)(c,v l )f 
Jy'i c 
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(49) 



T alS can also be rewritten as 
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I m 4 7 (w) 6 (c,^)(c,^')/rf 3 v. 
Jr 3 



(50) 



In eqs. (49) and (50), Y 3 is velocity space stretched by {y 3 \(-c, -c, -c) < (v 1 ^ 2 ^ 3 ) < (c, c, c)}. 
For convenience, non-dimensionalization is done as 
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(51) 



where L is the representative length in the observer's frame. 

With these non-dimensionalized quantities defined in eq. (51), the Maxwell-Jiittner function in eq. 



(30) can be non-dimensionalized as 



7(0) = ( mc ) 3 r(0) = n C -C7(")7(«)(i-S i C i ) 
1 noo 3 4vr^ 2 (C) 

(52) 

To solve eq. (47), the second-order TVD (Total Variable Diminishing) scheme |lUj is used for the left 
hand side of eq. (47) , and second order Runge-Kutta time integration is used for the time integration 
of eq. (47). 

The left hand side of eq. (47) represents the propagation of molecules with velocity vector (v 1 ,v 2 , t> 3 ) 
in physical space (x 1 , x 2 , x s ). This formulation for molecular propagation in physical space does not 
involve relativistic effects. As a result, in body-fitted curvilinear coordinates (£*, £ 2 , £ 3 ), molecules 
with a velocity vector (v 1 ,v 2 , t> 3 ) in (x 1 ,^ 2 ,^ 3 ) propagate with velocity Eq. (47) can 

therefore be written in body-fitted curvilinear coordinate (^ 1 ,^ 2 ,^ 3 ) as 



ll + ^d= (tWM*) + (c, -v'f x ^ 2 fj . (-vi) 

ot a£* \ ne + p ) c 1 t 

f = f/l 

where J is the Jacobian between x l and £\ 

The wall condition must also be considered. In this paper, complete diffusion at the wall is assumed. 
From conservation of the mass flux to the wall and by setting the £ 2 axis as the normal vector to the 
plane element of the wall, we obtain the following relation 

U = f(vl<0), (54) 
U/n w = f^(l,6 w ,0) (w|>0), (55) 

-Lf<o^f/r^ 3 v 

where f w is the distribution function on the wall. n w is the number density reflected from the wall 
and W is the temperature of the wall. 



V. RELATIVISTIC SHOCK LAYER PROBLEM 



In this section we consider the formation of shock layers around circular cylinders. We model this 
problem using the Anderson- Witting, Marie, and energy-preserved Marie model, which is obtained 
by setting 9 = 9e in eq. (38). For easier comprehension of physical conditions, for the observer's 
frame the absolute standard of rest is used as the hypothetical inertial frame. We use (x, y, z) instead 



of (x 1 ,^ 2 ,^ 3 ) and (v x ,v y ,v z ) instead of (v 1 ,v 2 ,v 3 ). The velocity corresponding to uniform flow is 
u x = 0.5c, u y = 0, u z = 0. Figure 2 shows a schematic view of the observer's frame and the flow field. 
The temperature of the uniform flow is #oo = mc 2 /45/c. Under these conditions, the Mach number of 
uniform flow calculated from eq. (33) is 2.689. The temperature of the wall is 9 W = mc 2 /30k. From 
eq. (32) r = A-Ka/L 2 = 10. Molecules are all assumed to be monatomic hard-sphere molecules. For 
the numerical grid, (v x , v y , v z , x, y) = (64,64,64,81,60). Numerical tests indicate that this numerical 
grid provides accurate simulations. Figure 3 shows the number density profile, the velocity profile, 
and the temperature profile along the stagnation streamline. The shock thickness simulated with the 
Anderson- Witting model is thinner than that simulated with either the Marie or the energy-preserved 
Marie models. The shock thickness simulated with the Marie model is thinner than that simulated 
with the energy-preserved Marie model. The difference of the simulated shock thickness between 
Anderson- Witting and Marie models is discussed later in detail by considering the behavior of the 
negative velocity tail of the distribution function. The temperature simulated with the Marie model 
is lower than that simulated with the energy-preserved Marie model. This implies that the dynamic 
pressure w in eq. (46) is negative and that energy dissipation occurs. 

Figure 4 shows profiles of heat flux q° and q x along the stagnation streamline. Both q° and q x have 
minimum values at the same point in the shock structure in all models and satisfy q a U a = in eq. 
(19). With the relation q a U a = 0, q° ~ near the boundary layer is demonstrated by the fact that 
the flow velocity is about zero for the stagnation point. q x is approximated from eq. (28) by using the 
N.S.F law. From the spatial gradient of 9 and p, qffSF approximated by the N.S.F. law is introduced 
from eq. (28). We define the heat flux by the gradient of the temperature as qf^gp = XV a 9 and by 
the gradient of the isotropic pressure as q X j^ SF = —X^j^Vp. As a result, q NSF = qffSF + Qn'sf- ^or 
results from the Anderson- Witting model, Figure 5 shows qfisF together with q x along the stagnation 
streamline on the left y axis, and on the right y axis Fig. 5 shows \Qn'sf/ ( 1nsf\ a l° n g t ne stagnation 
streamline, which is the ratio of absolute values of qfjsF anc ^ Qn°SF' ^ ne t nerma l conductivity A from 
the Anderson- Witting model, which is necessary for the calculation of q%sFi Insf an d Qn'sF' 1S & ven 
by pp. As shown in Fig. 5, q x < q%sp near 1.8 < —X/R < 3.2 in the shock structure indicates 
significant effects by the terms from the Burnett equation. Near —X/R ~ 1.4, Iq^sf/QnsfI exhibits a 
maximum value. As shown in Fig. 5, the heat flux calculated by the gradient of the isotropic pressure, 
Qn'sF' 1S nonnegligible for the calculated heat flux q x in this problem. 

Figure 6 shows the profile of the dynamic pressure along the stagnation streamline. Under steady 
flow conditions, w = u% Jj^ is obtained from eqs. (20) and (26). This indicates that the dynamic 
pressure depends on the product of the gradient of the number density and the flow velocity. As shown 



in Fig. 6, the dynamic pressure is negative in the shock structure and in the boundary layer, where both 
the gradient of the number density and flow velocity are positive. The dynamic pressure is positive 
in the sandwiched region defined by the shock structure and the boundary layer, 1.2 < —X/R < 2.0. 
Assuming that the approximation by the N.S.F. law in eq. (26) is adequate to describe the behavior 
of the dynamic pressure in both the shock structure and in the boundary layer, eq. (26) indicates 
that T) is negative for all models in the shock structure and the boundary layer. We are continuing 
to investigate whether or not this negativity is caused by contributions from the Burnett terms [9]. 
We do know, however, that the difference between the Marie and the energy-preserved Marie models 
in the bulk viscosity shown in Fig. 1 is not reflected in the difference of the profiles between these 
models, as shown in Fig. 6. However it is notable that the dynamic pressure has negative profiles 
in the boundary layer, where it might be less affected by the Burnett terms than in the shock structure. 



Figure 7 shows simulated distribution functions at selected points on the stagnation streamline 
for the Anderson- Witting and Marie models. In the shock structure, nonequilibrium conditions exist 
in the negative velocity tail, as shown in Figs. 7A and 7B. Near the positive peak of the dynamic 
pressure —X/R = 1.978, Fig. 7C indicates that nonequilibrium conditions exist near the peak of the 
distribution function. Fig. 7D indicates that at —X/R = 1.246, which is ahead of the boundary 
layer, that the distribution function represents a weak nonequilibrium condition at its peak. Fig. 7E 
indicates that at —X/R = 1.082, which is the middle of the boundary layer, that the distribution 
function shifts slightly to the left from the equilibrium distribution function. Fig. 7F indicates that 
at —X/R = 1.0, that the distribution function is not contiguous on both sides of v x /c = 0. 



Finally, we consider the relaxation process of the distribution function in the shock structure for 
the Anderson- Witting and Marie models by focusing on the dynamics of the negative velocity tail of 
the distribution function. To explain the numerical results shown in Figs. 7A and 7B, the relaxation 
rate for both models is clarified. 



We introduce a new relaxation rate parameter, 4>aw, which is obtained by removing n ^, +p and r 

(7(u)(l,«i)+ 1 &=)(l -v i ) T 

from the relaxation rate " e PJ of the right hand side of eq. (47) as 



^ AW = n — (57) 

(1 — u z ) ' 

We restrict ourselves to 4>aw along the stagnation streamline, for which Uj = (u x ,0, 0) and V{ = 



(v x ,v y ,v z ) gives 4>aw as the specific case of 4>aw, 

1 - u x v x 

* AW = V^F' m 

For comparison, the relaxation rate parameter for the Marie model, 4>m, is obtained by rewriting eq. 
(38) into H + = (/ (0) - /) and replacing 4> M = with v y = v z = from the right-hand 
side of eq. (47). Finally 4>m can be expressed as 

0M = (1-^) 1/2 >(1-^ 2 ). (59) 

Figure 8 shows 4>aw versus v x for various values of u x , and also shows (pM versus v x . When the gas is 
at rest, u x = 0, 4>aw = 1 f° r all values of v x . Generally 4>aw is the envelope for the half circle of 4>m, 
which contacts the half circle of 4>m at v x = u x . Because the absolute value of u x is approximately 
equal to the speed of light, 4>aw is approximately infinity, except for when \u x \ = c, in which case 4>aw 
approaches zero as \u x \ — > c. On the other hand, 4>m decreases as the absolute value of v x increases, 
regardless of the magnitude of u x . In the shock layer, the extent of nonequilibrium in the negative 
velocity tail determines the shock thickness, because nonequilibrium first appears in the negative ve- 
locity tail. Figure 7 indicates more rapid relaxation in the negative tail in the Anderson- Witting model 
than in the Marie model. The more rapid relaxation yields lower population Lower population of what? 
in the negative velocity tail. Consequently the shock thickness simulated with the Anderson- Witting 
model is thinner than that simulated with the Marie model. 



VI. CONCLUSIONS 



In this paper, we used two different relativistic Boltzmann-kinetic equations, the Anderson- Witting 
and Marie models, to numerically solve a shock-layer problem. The simulated heat flux has similar 
characteristics to that approximated by using the N.S.F. law. The heat flux calculated by using the 
gradient of the isotropic pressure represents a nonnegligible component of the total heat flux calculated 
for this problem. On the other hand, the simulated behavior of the dynamic pressure shows opposite 
tendencies to that approximated by using the N.S.F. law. It is important to determine whether or 
not this discrepancy between the N.S.F law and the simulated dynamic pressure is caused by the 
contribution of terms above the Burnett equation. The relaxation rate of the distribution function by 
using the Anderson- Witting model depends on both the flow velocity and on the molecular velocity. 
On the other hand, the relaxation rate simulated with the Marie model depends only on the molecular 
velocity. This difference of relaxation rate between these two models is manifested as a difference of 
the thickness of the shock structure. 



APPENDIX 



V>(Ce) in eq. (41) is given by 

^(&) = rmfcc 2 2 ° Gi; + 3fe " 13G|C f ~ 2 55 ! = + 2G|Cl , (A.l) 

(,eC v (C,e) 

where C v (Qe) is the constant-volume specific heat given by C v (Ce) = k (C, E + 5Ge(e — G%Ce ~ 
The correct bulk viscosity for 6 7^ Be is given by pQ 

, r PE k 2 {20G E + 3(e - 13G|Ce - 2G E Q% + 2G%Q E ) (4 - (% - 5G e (e + G E ( 2 E ) , A , 

" (Cb) = -3 gaceY (A - 2) 

The bulk viscosity for 6 = 6e is given by 

_ rpk 2 (20G + 3C - 13G 2 C - 2GC 2 + 2G 3 C 2 ) (l - C 2 - 5GC + G 2 ( 2 ) 
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FIGURE CAPTIONS 



FIG. 1: Bulk viscosity for the Marie model for 9 = 9e and 9 ^ 9e- 
FIG. 2: Schematic view of observer's frame and flow field. 

FIG. 3: Number density, temperature, and velocity profiles along the stagnation streamline for the 
A.W. (Anderson- Witting), Marie, and energy-preserved Marie models (9 = 9e)- 
FIG. 4: Heat flux (q° and q x ) profiles along the stagnation streamline for the A.W. (Anderson- Witting), 
Marie, and energy-preserved Marie models (9 = 9e)- 

FIG. 5: Heat flux q% SF approximated by using the N.S.F. law and Iq^sf/QnsfI along the stagnation 
streamline for the A.W. (Anderson- Witting) model. 

FIG. 6: Dynamic pressure profiles along the stagnation streamline for the A.W. (Anderson- Witting), 
Marie, and energy-preserved Marie models {9 = 9e)- 

FIG. 7 (A)-(E): Distribution functions and equilibrium distribution functions at selected points (A) 
-X/R = 4.965, (B) -Xj R = 3.218, (C) -X/R = 1.978, (D) -X/R = 1.246, (E) -X/R = 1.082, and 
(F) —X/R = 1.00 (Wall) on the stagnation streamline for the A.W. (Anderson- Witting) and Marie 
models. 

FIG. 8: Relaxation rate parameters (pM and 4>AW f° r various values of u x /c versus v x /c. 




FIG. 1 Bulk viscosity for the Marie model at B — Be and 8 ^ Be- 
(a is the collision cross section for hard-sphere molecules) 




FIG. 2 Schematic view of observer's frame and flow field. 



Number density, Velocity and Temperature profiles along the stagnation streamline 




FIG. 3 Number density, temperature, and velocity profiles along the stagnation streamline. 



Heat flux profiles along the stagnation streamline 
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FIG. 4 Heat flux (q° and q x ) profiles along the stagnation streamline. 



Profile of the heat flux approximated by N.S.F. law and lq X N SF /q x N s F ' a l° n g tne stagnation streamline 
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FIG. 5 Heat flux qff SF approximated by using the N.S.F. law and Wn'sf/InsfI a l° n g the stagnation 

streamline for the A.W. model. 



Dynamic pressure profiles along the stagnation streamline 




Distribution functions at -X/R=-3.218 
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FIG. 7 (B): -X/R = 3.218 



Distribution functions at -X/R=- 1.978 
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FIG. 7 (C): -X/R= 1.978 



Distribution functions at -X/R=- 1.082 




Distribution functions at -X/R=-1.000(WALL) 




Relaxational rate parameter vs v x /c for A. W. model and Marie model 
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FIG. 8 Relaxation rate parameters </>m and 4>aw for various values of u x /c versus v x /c. 



